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Abstract 

We study an extension of the supersymmetric standard model with four families and 
gauged horizontal symmetry SU(4)h, in which R parity automatically follows as a con- 
sequence of gauge invariance and of the field content of the model. Beyond ensuring R 
parity conservation, the horizontal symmetry allows for the dynamical generation of a 
hierarchical pattern of fermion masses, without the need of any ad hoc choice of small 
Yukawa couplings. The structure of the mass matrices implies that the fourth family does 
not mix with the other three, and that the b f , t', r', v' masses are all naturally in the 100 
GeV range. The model has other interesting phenomenological implications. The scale 
of the horizontal symmetry breaking is constrained by cosmological and astrophysical 
arguments to be ~ 10 11 GeV. Then the late b' decays could explain a magnitude of the 
isotropic cosmic gamma flux. In addition, a lower bound of a few eV is found for the mass 
of the r-neutrino, which can thus provide a hot component for the cosmological dark mat- 
ter. Due to R parity conservation, the lightest supersymmetric particle is stable and can 
provide the cold component. The v e and masses can naturally be in the correct ranges 
for a solution of the solar or atmospheric neutrino problems via neutrino oscillations. 
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1. Introduction 

In the Standard Model (SM), Baryon (B) and Lepton (L) numbers are conserved as a 
result of accidental global U(1)b and U(1)l symmetries that follow from the requirement 
of gauge invariance and renormalizability.Q In the supersymmetric standard model (SSM) 
this is not true anymore. In the superpotential, the renormalizable terms 

lle c , lqd c , u c d c d c . (1) 

(/, e c and q, u c , d c are the lepton and quark superfields) are allowed by the gauge symmetry, 
which violate L and B. In particular, the combination of first and third of these terms 
would lead to catastrophically fast proton decay mediated by <i°-type squark exchange. 
The relevant symmetry that ensures the absence of the terms ([!]) is called R parity, which 
is defined as R = (— i) 2J + 3B + L ) where J is the spin of the particle and B(L) its baryon 
(lepton) number f2|. R parity does not commute with supersymmetry (SUSY). On the 
other hand, it is an automatic consequence of a Z 2 matter parity under which the fermion 
superfields change the sign while the 'Higgs' ones $1^ remain invariant. 

An unsatisfactory feature of the SSM is that the Z 2 (or equivalently R) parity conser- 
vation has to be imposed by hand. The SUSY SU(5) model does not differ much from 
the SSM, since the gauge invariant term 10 5 5, involving the fermion superfields in the 
10 + 5 representation of SU(5) leads again to the set of B and L violating couplings ([!]). 
the SO(10) GUT offers an elegant solution to this problem: the Higgs fields are generally 
assigned to vector representations (10, 45, 54, 126 etc.) while the fermion superfields are 
in the spinor representation 16. While all the needed mass terms are allowed, the terms in 
(|I|) would arise from 16 3 which is forbidden by the gauge symmetry. In other words, since 
50(10) invariance allows for only pairs of 16-plets, the theory has an automatic Z 2 matter 
parity under which 16-plets change the sign whereas the superfields in vector representa- 
tions remain invariant. However, this is not true anymore for the 5*0(10) models in which 
symmetry breaking is triggered also by the scalar components of superfields belonging to 
the 16$ + 16$ [||]. In fact, after substituting the VEV (16$), the couplings -pl6 3 16$ 
lead again to the terms ([l]). Since in these models the right handed neutrino masses are 
generated by operators ~ -gl6 2 16| of the same structure, the ratio (16$)/M cannot be 
too small, implying in turn that the magnitude of the resulting R-parity violating terms 
conflicts again with the limits on the proton lifetime. Thus R-parity conservation is not 
automatic in 50(10) models with Higgses belonging to the 16$, and matter parity has 
to be imposed by hand in order to distinguish the fermion 16-plets from the Higgs ones. 

In this paper we wish to put forward the idea that R (or equivalently Z 2 ) parity 
conservation can be naturally ensured in models based on gauged non-Abelian horizontal 
symmetries. Such models are particularly interesting since they can also explain at least 
qualitatively the observed pattern of fermion masses and mixing. Namely, the structure 
of the fermion mass matrices can be related to the horizontal symmetry breaking pattern, 
and the mass hierarchy between families arises dynamically from certain hierarchies in this 

1 These symmetries can be broken only by non-renormalizable Planck scale operators like jjUQQ Q. 
This term violates L and induces the neutrino Majorana masses of about 10 -5 eV, which could be relevant 
for the solar neutrino oscillations. However, analogous terms violating baryon number are very small to 
cause any observable effect. 
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breaking (see for example the models P, [5] based on SU(3)h symmetry). In this kind of 
models it is natural to assume that the horizontal group Gh acts only on the quark-lepton 
superfields, while the Higgses $1,2 responsible for the electroweak symmetry breaking are 
Gtf-singlets. Hence, independently of the choice of the vertical gauge group and/or of the 
particular superfield assignments to its representations, the Higgs and fermion superfields 
can be always distinguished since they carry different horizontal quantum numbers. This 
leads to the possibility of allowing the necessary mass terms which are bilinear in the 
fermion superfields, while forbidding the trilinear couplings ([!]). 

Our task is now to find and classify the theories in which the horizontal gauge group 
Gh automatically forbids the unwanted B and L violating terms in ([!]). Then R parity (or 
equivalently Z 2 matter parity) will appear as an automatic consequence of the horizontal 
gauge symmetry and of the field content of the model. We demand that the models we 
are interested in should satisfy the following list of basic requirements: 

(i) In order to ensure a straightforward definition of the horizontal gauge symmetry, 
all the fermion superfields q, I, u c , d c and e c should be assigned to the same representation 
of the non-Abelian horizontal group Gh- In other words, we forbid Gh singlet families. 

(ii) We require that the couplings in ([D are forbidden by gauge invariance. Therefore 
SO(3)h and chiral SU(3)h, which are the only simple groups containing three dimensional 
representations, are excluded since in both cases 3 3 contains a gauge singlet. However 
we immediately notice that for SU(N)h (N > 3) the term N 3 does not contain gauge 
singlets. Hence these groups represent a class of interesting candidates. 

To have phenomenologically realistic theories, the following additional constraints 
should be also imposed: 

(Hi) In order to avoid the proliferation of Higgs doublets with masses at the electroweak 
scale, the standard Higgs superfields $1^ must be Gh singlets. The presence of several 
Higgs doublets would in fact spoil the natural suppression of Flavor Changing Neutral 
Currents (FCNC). It would also destroy gauge coupling unification, thus preventing any 
attempt to embed the model in some vertical GUTs. 

(iv ) The model should provide naturally a realistic pattern of the fermion mass and 
mixing. In this respect, models based on chiral SU(3)h @, U have proven to be quite 
successful in relating the fermion mass hierarchy to a hierarchy in the horizontal symmetry 
breaking VEVs, and in accounting for the observed fermion mass and mixing pattern. 
Even if chiral SU(3)h fails to satisfy condition (ii), this again suggests that chiral SU(N)h 
(N > 3) are promising candidates.^ 

(v ) A final strong condition is that R-parity breaking terms should not appear even 
after Gh breaking. In other words, no effective couplings which could generate the terms 
in ([D after substituting the VEVs £ fe — > should be allowed by the Gh symmetry. 
In particular, this condition restrict the viable SU(N)h models to even N. Consider 

2 We note that vectorlike SU(3)h with q, I transforming as 3 and u c , d c , e c as 3, forbids the first two 
terms in d), which is enough to ensure proton stability. However it allows for the SU (3)h invariant 
Yukawa terms // c $i,2- Then the mass splitting between different families can be achieved only at the 
price of several unnatural fine tunings between horizontal non-invariant effective operators. Thus this 
case can be hardly regarded as realistic. (In addition unification of the quarks and leptons within one 
irreducible GUT multiplet would clearly be not possible). The same problem is encountered for SO (A) 
and 50(5) with 4 dimensional representations. While the 4 3 terms are forbidden and condition (ii) is 
fulfilled, the invariant Yukawa terms // c $i,2 are again allowed, and thus condition (iv) is not satisfied. 
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in fact SU(N)h with the / and f c fermion superfields assigned to the fundamental N 
dimensional representation. The mass terms transform as N x N and thus belong to 
two-index (symmetric and antisymmetric) representations. In order to construct horizon- 
tal gauge invariant mass terms we can take also the horizontal Higgses £ k in two-index 
representations. Then for iV=4,6,. . .no invariants of the form N x N x N x V(£) (with 
V some polinomials in the horizontal fields £ fc ) are allowed. In contrast, for N odd the 
totally antisymmetric e tensor allows to rewrite some combinations of Higgs fields with 
an even number of free indices as tensors with an odd number of free indices, suitable 
for generating effective N x N x N x V(£) gauge invariant terms. After the breaking 
of the horizontal symmetry (£ fc — > (£ fc )) these effective operators will then produce the 
R- violating terms in (|l]). 

This brief analysis suggests that natural conservation of R-parity could be achieved 
in models based on chiral horizontal symmetries SU(N)h with N even, under which 
the quark and lepton superfields transform as fundamental iV-plets. Hence the number 
of families must be extended to Nf = 4,6,.... As is well known, the possibility of 
extra families with a light neutrino is ruled out by the measured invisible decay width 
of Z-boson, however, sequential generations with heavy neutrinos (m > Mz/2) are not 
excluded. On the other hand, detailed studies |J of the effects of radiative corrections 
due to additional families show that precise electroweak data are compatible with a fourth 
family, while six families (which would be the next interesting case) are ruled out. In 
addition, a dedicated analysis showing the viability of SUSY models with four families 
with respect to gauge coupling unification was presented in ref. f7j. These results are 
relevant for our analysis, since condition (Hi) ensures that the field content in our SU(4)h 
model is the same than that of the four family SSM of ref. , up to some large energy 
scale where the horizontal symmetry breaks down (see Sect. 2). Hence, we conclude 
that if natural R parity conservation arises due to some non-Abelian horizontal gauge 
symmetry, then the theoretical and phenomenological constraints hint to models based 
on the SU{4)h group, on which we will concentrate in the rest of the paper. 



2. SU(4)h symmetry and its phenomenological implications 

Let us consider the standard SU(3)xSU(2)xU(l) vertical gauge group, with local 
chiral SU(4)h horizontal symmetry acting on four families of left chiral superfields 

/«: ? * = (d) ~( 3 >2,l/6)«, L= (l) "O-' 2 '- 1 / 2 ^ 

fa- < ~(3,l, Q -2/3) Q , t£~(3,l,l/3)«, %~(l,l,l) a (2) 

where each superfield is assigned to the fundamental 4 representation (a = 1, . . .4 is the 
SU(4)h index). With this field content the horizontal SU(4)h is anomalous. In order to 
cancel the horizontal anomaly we introduce the following superfields which are vectorlike 
with respect to SU(3) x SU(2) x U(l) and belong to the 4 of SU(A) H : 

F a : U a ~ (3,1,2/3)", D a ~(3 i l,-l/3) a , E a ~ (1, 1, -l) a 

F c a : U?~ (3, 1,-2/3)", ^-(3,1,1/3)", (1,1,1)°, JV?~ (1,1,0)° (3) 
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As we will see in short, these superfields turn out to be necessary also for providing masses 
to the known fermions. 

In the Higgs sector, we choose the standard Higgs doublet superfields $1,2(1, 2, =Fl/2) 
to be singlets under SU(4)h- In order to break the horizontal symmetry and to generate 
the fermion masses, we introduce also a set of 'horizontal' Higgs superfields £{ Q/ 3} (k = 
1,2 . . .) transforming as the symmetric 10 representations of SU(A)h (clearly, they must 
be SU(3) x SU{2) x U(l) singlets). Additional superfields £ fc transforming as 10 are also 
needed to render the Higgsino sector free of chiral anomalies. With this field content, the 
most general Yukawa terms allowed by gauge invariance read as 



W F = g f f a F c a $ 1(2) + £ h k F F a F? Hp + kfF a f c a , F = E,D, (U) (4) 

k 

W N = gJ a N? $ 2 + £ hfrNZN! &p (5) 

k 

where the g's and Ks are Yukawa couplings which we take to be all 0(1). The last term 
in eq. (f|) is a gauge invariant bilinear, and the Aj's are some large mass parameters. 
As already stated, no terms trilinear in the quark and lepton superfields are allowed by 
the SU(4)h gauge symmetry, ensuring naturally the absence of the B and L violating 
couplings (U). We are facing here a situation analogous to the 50(10) model, since Im- 
parity does not have to be imposed by hand, but appears as an accidental symmetry that 
follows from horizontal gauge invariance. The terms (||) are in fact invariant with 
respect to the Z 2 parity under which the fermion superfields /, f c , F, F c change sign, while 
the Higgs superfields $1,2 and £ stay invariant. More in general, these have an automatic 
global symmetry U(1)h under the following transformations: 

fj c ^e^f,f c , F,F c ^e~^ F, F c , - e 2 *^*, $ 12 -> <J> 1)2 (£* -> e^ k ) (6) 

Its Z 2 subgroup (u = it) remains unbroken even when the scalars £ get non-zero VEVs. 
Clearly, this Z 2 matter parity which ensures R parity conservation and hence proton 
stability will be respected by all possible higher order Planck scale operators as well. 

The Yukawa terms (f|) lead to the so called "universal seesaw" mechanism 0, |9| for 
generating the masses of the charged fermions, while the terms (|5]) lead to the ordinary 
seesaw mechanism for neutrino masses ||lQj|. After the horizontal scalars £ fc develop non- 



zero VEVs, the extra fermions F and F c of eq. (j^) acquire large masses through the 
second term in eq. (ffl). Then the first and third terms cause a "seesaw" mixing of the 
ordinary quarks and leptons /, f c with the heavy ones. As a result, in the base (/, F) 
(f c , F c ), the 8x8 mass matrix for the charged fermions f — e,d, (u), F = E,D, (U) reads 

m, = (? 9, y), < = E^O (T) 

k 



M F 



where V\i = ($1,2) are the VEVs of the two electroweak Higgs doublets. As for the 
neutral states, in the base (u, N c ) the 8x8 Majorana mass matrix has the form 

Mv = {L P)> ^5 = E^<&>- (8) 
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In contrast to the SM and most GUT models, in our picture the fermion mass hierarchy is 
not generated by an ad hoc choice of the Yukawa couplings. In fact we can assume all the 
Yukawa constants to be 0(1), close to the size of the gauge couplings. As long as the off- 
diagonal blocks in eqs. fl7|) and are flavour blind (unit) matrices, all the informations 
on the fermion mass and mixing pattern is contained in the heavy mass matrices M F . 
The structure of the latter is determined by the different VEVs (£ fc ) (modulo differences 
in the Yukawa constants hp), and the observed hierarchy of the light fermion masses is 
ultimately determined by the hierarchy in the VEVs which break the SU(4)h symmetry. 
In other words, the VEV pattern should provide a step-by-step breaking of the chiral 
horizontal symmetry 

SU{A) H x 1/(1)* 3 SU(3) H x U(1)' H 3 SU(2) H x U(1)» H 3 £7(1). (9) 

The first breaking at the scale V\ (given by the VEVs (£11)) provides the mass terms for 
the first heavy family F\, the second breaking (at V 2 ~ (£^ 2 )j a — 2) f° r the second family 
F 2 and the third stage of breaking at V 3 ~ (6^3) > a < 3 generates the mass of the F 3 
state. At this stage all the horizontal gauge bosons have acquired large masses, the less 
massive are the ones responsible for the flavour- changing transitions between the fourth 
and the third families, with masses ~ V3. The model also provides a natural possibility 
to obtain together with three light families, a heavy fourth, for which all the masses are 
of the order of the electroweak scale. In fact let us assume that all the VEVs of the type 
(^4) are vanishing, so that a diagonal global £7(1) subgroup of SU(4) H x U(1)h, given by 
the generator T = diag(0, 0, 0, 1), is left unbroken.^] Then the 4x4 mass matrices M F ( N ^ 
are rank-3 matrices of the following form 

^ F= ( M i 0) ' f = u > d > e > n ( iq ) 

where the 3x3 blocks M F contain non-zero entries, given by the VEVs of the horizontal 
scalars. Clearly, the residual global symmetry £7(1) ensures naturally the heaviness of the 
fourth family, since there is no seesaw mechanism for the corresponding fermions. The 
right-handed components of the fields f'4 = b', t', r' , v' are actually the states, whereas 
the fl form with the F4 superheavy particles of mass A/. From eqs. (0) and (§) we obtain 

my, T / = gd,eVi, rrify = 9u,yV 2 - (11) 

Since the Yukawa couplings are all 0(1), for moderate values of tan/3 = v^/vi all the 
masses in ( [TT]) are of the order of electroweak scale ~ 100 GeV. On the other hand, 



it is apparent from (0) that the fermions of the first three families will acquire their 
masses through a seesaw mixing with the superheavy F fermions. Their masses will 
then be suppressed from the electroweak scale down to the observed values, provided that 
Mp^ 1 > Af. Namely, after decoupling the heavy states^ the 3x3 mass matrices of the light 

3 A detailed analysis of the horizontal Higgs potential will be given elsewhere. 

4 Below the scale Vh — V3 our theory simply reduces to the SSM with four families. Therefore, for 
Vh large enough all FCNC phenomena related to the horizontal gauge or Higgs bosons and to the mixing 
with weak isosinglet heavy fermions are strongly suppressed. 
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charged fermions / = e,d, (u) are ml 3 ' = gf/J,f(Mp) 1 fi(2) while for the 3x3 Majorana 

mass matrix for the light neutrinos we obtain = (M^ )~ l (g v V2) 2 ■ Hence through the 
seesaw mechanism the horizontal VEV pattern (H) is reflected in the observed pattern of 
fermion masses in an inverted way || (see also |JI3"]). Modulo the different 0(1) Yukawa 
constants hp the heavy mass eigenstates of M F reflect the hierarchy V\ ^> Vi V3 
of the SU(4)h symmetry breaking. On the other hand, the ordinary quark and lepton 
masses are inversely proportional to the masses of their heavy partners. For example, for 
the charged lepton and the neutrino masses we haveQ 

^ A e m T / ^ ml, 

m e,fi,T — 1^3 < m r' 1 m f e ,u l j„u T — ,1^3.. ^ m u' (t^J 

hp ^1,2,3 If'N ^1,2,3 

The quark masses are given by expressions analogous to that of charged leptons. The CKM 
matrix of mixing between the first three families is also determined by the structure of the 
matrices Mp\ and arises from the presence of non-diagonal VEVs (££g) (en, ft = 1,2,3). 
The fact that the b and r masses are of order a few GeV, implies that the masses of the 
corresponding heavy states h z D V^ and hpV^ are not much larger (say, within one or two 
orders of magnitude) than the mass scale A^ e . As for the top quark, the value of its mass 
m t ^ 150 GeV requires hfjVs ~ A u . In this case corrections to the seesaw formula should 



be taken into account in relating mt to the heavy scales (see e.g. ||13|| ). 

As we stated above, the residual U(l) symmetry implies that fourth neutrino v' is a 
Dirac particle with mass ~ 100 GeV, while three neutrinos get small Majorana masses 
given by eq. flT2|). Therefore, the neutrino mass hierarchy is expected to be approximately 
the same as the hierarchy between the quarks or the charged leptons: 

m Ue : : m Ur ~ m u : m c : m t or ~ m e : m M : m T . (13) 

On the experimental side, the firmest constraints on the masses of any new sequential 
fermion, quark or lepton, have been set at LEP: m/ ^ Mz/2. This indeed represents 
the best constraint on m T i and m u > . Searches for new quarks at the TEVATRON collider 
could in principle give much better bounds for m t i and my [PHI . However, let us note that 
the structure of the heavy mass matrix ( |10"1) implies that the fourth family is unmixed 
with the three lighter ones. Hence the usual signatures, as for example b' —>■ c,u, that 
have been used to set the limits on new sequential quarks [|TTJ do not occur in our case. 



In the absence of a detailed experimental analysis of the unmixed case, the only reliable 
limit is again the LEP one also for the new quarks. Hence we can safely conclude that 
the predictions in ( |TT| ) are by no means in conflict with the existing experimental limits. 
However, it is clear that for the masses of the fourth family fermions not much room is 
left. The allowed parameter space is strongly constrained by the CDF measurement of 
the top mass, mt = 174 ± 10 ± 13 GeV |12| , by the precision tests of the SM which do not 
leave much space for additional sizeable radiative corrections as would be induced by a too 
large m t '-my splitting, and by renormalization group analysis of the Yukawa couplings, 
much in the spirit of ref. [[/]]. In particular, the universal seesaw mechanism implies 
m t i > m t , and most likely m t > > my. Then, according to H, for m t i > m t > 150 GeV the 

5 Since in the following we deal with order of magnitude estimates, we neglect the different 0(1) factors 
accounting for the renormalization group running of masses from the horizontal scale to lower energies. 
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consistency of the model implies not too large masses for the fourth family fermions. For 
the low values of tan/3 we are interested in (e.g. tan/? ~ 2), the maximal values allowed 
are about my ~ 100 GeV and m T /y ~ 50 GeV, that is within the reach of LEP II. 



3. Astrophysical consequences of the model 

Indeed the presence of a fourth family of fermions constitutes the most compelling 
prediction of our model, so we will now address some issues regarding these new states. 
The unbroken diagonal U{1) subgroup of SU(A) H x U{1)h implies that the fourth fam- 
ily is unmixed with the three lighter ones. We assume that the lightest member of the 
fourth generation is the neutral one z/, as is also suggested by the analysis of ref. 0. 
For simplicity we also assume that m t > > my. Then V and v' are stable with respect to 
electroweak interactions. A stable v' with mass in the 100 GeV range is cosmologically 
safe, the contribution p u > to the present cosmological density would not exceed the crit- 
ical density p c even in the presence of a v'-v' primordial asymmetry of the order of the 
baryon asymmetry. However, as we will argue in the following, in our model no sizeable 
asymmetry has to be expected for the fourth family fermions. 

In contrast, the existence of stable heavy quarks carrying colour and electric charge 
would constitute serious problem, since it will conflict with the constraints arising from 
superheavy element searches |15|]. Indeed, the stable b' would behave essentially as d 
quarks, hadronising into heavy stable mesons b'u and 'nucleons' b'ud [IIS] giving rise to 
heavy hydrogen-like 'isotopes' with masses ~ 100 GeV. The existing experimental limits 
on this kind of isotopes are extremely tight. For my < 1 TeV the limit on the V abundance 
relative to normal hydrogen is ny /tib < 10 -28 [16]. 

However, the exchange of the SU(4)h gauge bosons Zh would allow the heavy quark 
to decay, dominantly through the channel b' — > bu T u', with a lifetime 

where Vh = V3 is the lowest scale in the horizontal symmetry breaking (see eq. @), 
v = 174 GeV is the electroweak scale and r M = r(/i — > ev^^) is the muon lifetime. We can 
use cosmological arguments, together with the experimental limits on searches for heavy 
isotopes, to put an upper bound on Ty, which in turn will translate in an upper limit 
on Vh- Taking into account the finite lifetime of the heavy quarks, their present number 
abundance relative to baryons is ny/ns = ?"oexp(— t /ry) < 10~ 28 , where r = (ny/riB)o 
is the relic abundance for stable b', h = 0.5 — 1 is the Hubble parameter in units of 100 Km 
s _1 Mpc _1 , and to — 2 • lO 17 /!" 1 s is the present age of the Universe (we assume matter 
dominated expansion of the Universe, and Q = 1 as it is motivated by inflation). From 
this equation we get an upper limit on the b' lifetime 

ry < 3.1 • lO 15 /*- 1 ^ + 0.036 lgror 1 s. (15) 

Due to many theoretical uncertainties related to the actual annihilation cross section for 
the b', it is not possible to compute precisely the value of r 0r However, an estimate 



of the relic abundance of heavy stable <i-type quarks has been given in [15|. Under the 



7 



assumption that there is no cosmological b' and b' asymmetry, it was found that for 
my ~ 150 GeV the energy density of these relics, relative to critical density (namely 
Qyh 2 ) could range from 10 -9 to 10~ 4 (smaller values are obtained for lighter b' masses). 
The lower limit corresponds to the case when the relic density is determined by the 
annihilation after the QCD phase transition, and it was obtained by taking as an upper 
bound on the annihilation cross section the geometrical cross section (ctq ~ 100 mb ). 
The upper limit was obtained under the opposite assumption, namely that annihilation 
after confinement is negligible, and that the relic density is essentially determined by the 
QCD annihilation cross section for free quarks. Then the ratio of the b' to baryon number 
densities r = (ny /ub)o = fly /flB' m B/ m b' lies in the range 3- 10 -10 < r < 3- 10~ 5 where 
we have taken Q B ~ 0.02 as suggested by nucleosynthesis estimates. As is discussed in 



TB| , the most reasonable assumption is that the relevant annihilation process happens 
after confinement, however with a cross section much smaller than the geometrical one, 
giving ro ~ 10~ 7 — 10~ 8 . Clearly, in the presence of a sizeable baryon asymmetry in the 
fourth family sector, the relic abundance of the heavy V quarks would be some orders of 
magnitude larger than the quoted estimates, up to r ~ 1. 

As we see, the bound (]15[) very weakly depends on the initial b' abundance. Even if 
we let ro range between 1 — 10~ 10 , by taking h = 0.5 we obtain ry < 6 • 10 15 — 10 16 s. On 
the other hand, according to eq. (0), this bound translates into an upper limit 

Vu < /15 ° ■ 3 ■ 10" GeV < 4 ■ 1Q» GeV (16) 

- h l / A (l + 0.036 lgro) 1 ^ - v > 

where we have taken the conservative upper bound my < 150 GeV from the analysis 0. 
For more realistic values my ~ 100 GeV |7| we get Vh < 2.4 • 10 11 GeV. According to 
eq. fljjD , this upper limit on the scale Vh = V3 together with the experimental limit 
77V > Mz/2 translates into a lower bound on the r-neutrino mass: 

2 

m Vr ~ hjf 1 ^ £ (1 - 10) eV (17) 

'H 

where in the numerical estimate we have taken into account the 0(1) uncertainties in 
the Yukawa coupling (for perturbativity we have to assume h N < 3). On the other 
hand, the most conservative cosmological upper bound on the light stable neutrino masses 
m UT < 92 Q h 2 eV provides a lower bound on Vh' 



v " ^w?(^)^ (0 - 3 - r > ■ 1011 GeV < 18 > 

More stringent limits on ro and Ty can be derived by considering that the late decay 
of the b' can cause a significant contribution to observed cosmic ray fluxes, in particular 
to the isotropic diffuse gamma-ray background [17]. Indeed, at the moment of decay, 



the b' quarks are bounded within colorless hadrons like b'ud or b'u jl5fl . Then the decay 
b' — > bv T v' will produce an hadronic jet with the b quark being the leading particle and 
an excitation energy E ~ ^rriy. The fragmentation of this jet produces 7r°, r\ etc., with 
the subsequent radiative decay resulting in a specific photon spectrum, and the number 
of photons produced is directly proportional to r . In order to estimate the flux at the 
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present era, the redshift in the photon energies has to be taken into account as well. 
As long as the decay happens at the matter dominated epoch, and the small amount 
of relativistic decay products does not affect sensibly the Universe expansion rate, we 
have 1 + z — (to/ r b') 2 ^ 3 ~ 10 — 20 for the values of r&> estimated above. To compute 
the value of the isotropic cosmological gamma-flux we need to know what fraction of the 
jet energy E Q is taken by the photons and what is the energy spectrum. The photon 
spectra produced in jet hadronization for different leading particles were computed, using 
a Monte Carlo simulation program, in ref. [|18j. These spectra exhibit a remarkable 
scaling property in terms of the variable x = E^/E , and in the case the leading particle 
is a b quark, the photons carry away about 25 percent of the initial jet energy. Using 
the results of ref. [18| we have estimated the present value of the gamma-flux d^^/dE^ 



as a function of r , and we have compared it with the existing observational limits (see 
T9| and references therein). For example, for E 1 = 100 MeV the experimental bound 



is d<& 1 /dE 1 ^ 10~ 7 cm -2 s _1 sr _1 MeV -1 and the cosmic gamma-flux produced from b' 
decays at z = 10 — 20 saturates this bound for r ~ 10 _7 .f] Thus, for the preferred values 
ro ~ 10~ 7 — 10~ 8 [I5| our model is consistent with with the observational data on the 



isotropic gamma-background, and can even provide an appealing explanation of its origin. 

The previous analysis implies that r should be rather small. In particular this ex- 
cludes any sizeable baryon asymmetry between b' and V . Then the possible baryogenesis 
mechanisms applicable to our model are severely constrained. Even if a baryon asym- 
metry in the fourth family is hardly expected, since there is no mixing with the other 
three families, and hence no source of CP violation, the high rate of sphaleron processes 
before the electroweak phase transition |2I| would immediately redistribute any baryon 
asymmetry present in the first three families to the fourth one. Therefore, no mecha- 
nism is acceptable which generates the baryon asymmetry before the sphaleron effects are 
switched off.|] In the context of our model, the most appealing scenario is to assume that 
baryon asymmetry is generated at the electroweak (first order) phase transition, as a re- 
sult of CP violation in the dynamics of quarks or leptons interacting with the walls of the 
expanding bubbles of the broken phase p3f . Outside of the bubbles the electroweak sym- 



metry is unbroken, quarks are massless and the rate of the fermion number violation due 
to sphaleron transitions greatly exceeds the Universe expansion rate. Inside the bubbles 
the quarks are massive due to non-zero VEVs of the Higgs fields, the sphaleron processes 
are strongly suppressed and fermion number is effectively conserved. Baryon asymmetry 
inside the bubbles could be produced (and maintained) due to CP violating effects, as a 
difference between the quark and anti-quark fluxes penetrating the walls from the unbro- 



6 Substantially larger rg would require much larger redshift, and hence much smaller t^. However, 
the lower bound (|l8|) on the breaking scale Vh excludes much smaller lifetimes. 

7 In principle, in our model the baryogenesis with non-zero B— L could occur due to CP violation effects 
in out-of-equilibrium decays N c — > I + $ of the heavy right-handed neutrino [^lj (for the viability of this 
mechanism in the SUSY case see ref. 0|), or in the decays of SU(4)h gauge or scalar bosons. Then 
sphaleron effects would immediately transfer the produced net lepton number into a baryon asymmetry 
also in the fourth family sector. However, our model naturally avoids the possibility of such a lepto- 
baryogenesis. As it was shown in ref. |2^] , the large scale density fluctuations hinted by the COBE 
measurements require rather low inflationary reheat temperature (Tr ~ 10 s GeV) and correspondingly 
low inflaton mass (m n ~ 10 11 GeV). Then for Vh ~ 10 11 GeV the masses of all right-handed neutrinos 
and horizontal bosons are m„, and therefore they cannot be produced after inflation. 
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ken phase to the broken one. This will affect only the first three family fermions. Since 
the fourth family is unmixed and has no CP violation, no baryon excess is expected in 
this sector. Although the viability of such a baryogenesis in the SM is still disputed in the 
literature, in the context of SSM it could be more effective and sufficient for providing the 
observed baryon asymmetry. Clearly this topic deserves additional special considerations. 

Few remarks about the neutrino masses are now in order. For the r-neutrino mass 
our model predicts the lower bound eq. (|T7|) of a few eVs , and according to (|T3) 
and v e are expected to have much smaller masses. A v T with mass in the range 1 — 10 
eV will give a sizeable contribution to the cosmological energy density as a hot dark 
matter (HDM) component. We remind here that the COBE measurements of the cosmic 
microwave background anisotropy, together with other data on the density distribution 
of the Universe at all distance scales (galaxy-galaxy angular correlations, correlations of 
galactic clusters, etc.), can all be fit by assuming some HDM admixture to the dominant 
CDM component ||24|| . The best fits hint to a neutrino mass vtl Vt ~ 5 — 7 eV |25] which 
does appear naturally in our model. As for the CDM itself, in our R parity conserving 
SUSY model it can be naturally provided by the lightest supersymmetric particle (LSP), 
presumably a neutralino. 

As we commented earlier, the neutrino mass hierarchy should be qualitatively the 
same as that for the charged quarks and leptons. However, the spread in the Yukawa 
coupling constants hp does not allow to put severe limits on the other neutrino masses. 
For example, by taking m u ^/m UT ~ m c /m t , as is suggested by the first estimate in eq. 
(fjjl), one obtains m v ~ (2 — 5) • 10~ 3 eV. This range corresponds to the MSW solution 
of the solar neutrino problem |2S| via v e — > oscillations. Alternatively, if we had to 
attempt an explanation of the deficit of the atmospheric via z/ M — > v e oscillations (for 
a recent analysis, see |27|]), we would need m v ~ 0.1 eV which is compatible with the 
second estimate in eq. fll3|). Obviously the MSW explanation to the solar neutrino deficit 
would not be viable in this latter case. 



4. Conclusions 

In this paper we have put forward the idea that natural conservation of R parity in 
SUSY models can be guaranteed in the presence of a suitable horizontal gauge symmetry, 
so that an accidental Z 2 matter parity (equivalent to R parity) automatically follows from 
gauge invariance and the field content of the model. On theoretical and phenomenolog- 
ical grounds, we have identified SU(4)h as the only viable horizontal gauge group for 
implementing this idea. It implies existence of four fermion families. We have suggested 
a realistic SUSY model based on the SM vertical gauge group SU(3) x SU(2) x U(l) 
with an SU(4)h anomaly free horizontal gauge symmetry. It leads to a particular form 
of the fermion mass matrices rendering all the fourth family fermions naturally in the 
100 GeV range and unmixed with the first three families. This is due to the global U(l) 
symmetry which remains unbroken and in fact represents the conserved fermion number 
of the fourth family^] As for the masses of the light three families, our model leads to 

8 Let us remark that actually this is a rather general statement which applies not only to our particular 
SUSY model designed for automatic R parity: if the fourth family fermions will be indeed discovered in 
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"seesaw" suppression of their magnitude from the electroweak scale down to the observed 
values. This suppression is achieved dynamically, without the need of any tuning of the 
Yukawa couplings which have been assumed to be all 0(1). By means of cosmo logical 
and astrophysical arguments, we have constrained rather precisely the scale Vh at which 
the horizontal gauge symmetry is completely broken, obtaining a narrow window around 
10 11 GeV. The upper bound on the scale V H sets a lower limit on the r neutrino mass 
of about few eV. Hence our model naturally provides cosmological HDM in the form of 
v T and, due to R parity conservation, also CDM in the form of stable LSP. Since in our 
scheme conservation of R-parity is ensured by the horizontal gauge symmetry indepen- 
dently of the particular choice for the vertical gauge group, it would be interesting to 
extend the present analysis to phenomenologically appealing GUT models, such as SU(5) 
or E 6 , for which R-parity conservation is not automatic. 
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